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ABSTRACT 



In this work the double vector meson production in two-photon interactions at high energies 
is investigated considering saturation physics. We extend the color dipole picture for this 
process and study the energy and virtuality dependence of the forward differential cross section. 
Comparison with previous results is presented and the contribution of the different photon 
polarizations is estimated. 



1 Introduction 



The high energy limit of the perturbative QCD is characterized by a center-of-mass energy which 
is much larger than the hard scales present in the problem. In this regime the parton densities 
inside the projectiles grow as energy increases, leading to the rise of the cross sections. As long 
the energy is not too high, we have low values of the partonic density and the QCD dynamics 
is described by linear (BFKL/DGLAP) evolution equations [1, 2]. However, at higher energies 
the parton density increases and the scattering amplitude tends to the unitarity limit. Thus, 
a linear description breaks down and one enters the saturation regime, where the dynamics 
is described by a nonlinear evolution equation and the parton densities saturate [3, 4]. The 
transition line between the linear and nonlinear regimes is characterized by the saturation scale 
Qs&t{%), which is energy dependent and sets the critical transverse size for the unitarization of 
the cross sections. In other words, unitarity is restored by including nonlinear corrections in 
the evolution equations. Such effects are small for kj_ > Q 2 at and very strong for k\ < Q 2 at , 
leading to the saturation of the scattering amplitude, where k± is the typical hard scale present 
in the process. The successful description of all inclusive and diffractive deep inelastic data at 
the collider HERA by saturation models [5, 6, 7, 8] suggests that these effects might become 
important in the energy regime probed by current colliders. Furthermore, the saturation model 
was extended to two-photon interactions at high energies in Ref. [9], also providing a very good 



description of the data on the 77 total cross section, on the photon structure function at low 
x and on the 7*7* cross section. The formalism used in Ref. [9] is based on the dipole picture 
[10], with the 7*7* total cross sections being described by the interaction of two color dipoles, in 
which the virtual photons fluctuate into (For previous analysis using the dipole picture see, e.g., 
Refs. [11, 12]). The dipole-dipole cross section is modeled considering the saturation physics. 
The successful descriptions of the 77 interactions and light /heavy vector meson production for 
ep collisions at HERA are our main motivations to extend this formalism to describe the double 
meson production and analyze the effects of the saturation physics. 

In the last few years the double meson production has been studied considering different 
approaches and approximations for the QCD dynamics [12, 13, 14, 15, 16, 17, 18, 19, 20]. In 
particular, in our previous paper in Ref. [15], we have performed a phenomenological analysis 
for the double J/^f production using the forward LLA BFKL solution. In that case, the hard 
scale was set by the charm quark mass. There, we also studied the possible effects of corrections 
at next to leading approximation (NLA) level to the BFKL kernel investigating the influence of 
a smaller effective hard Pomeron intercept. Afterwards, in Ref. [17] the non- forward solution 
was considered for a larger set of possible vector meson pairs, where the large t values provide 
the perturbative scale. Moreover, in that paper the double vector meson production in real 
photon interactions was studied, the t-dependence of the differential cross section was analyzed 
in detail and the total cross section for different combinations of vector mesons was calculated 
using the leading order impact factors and BFKL amplitude. More recently, two other studies 
on the process 7*7* — > VV have appeared in literature [19, 20]. In the first one [19], the leading 
order BFKL amplitude for the exclusive diffractive two-p production in the forward direction is 
computed and the NLA corrections are estimated using a specific resummation of higher order 
effects. In the last paper [20], the amplitude for the forward electroproduction of two light 
vector mesons in NLA is computed. In particular, the NLA amplitude is constructed by the 
convolution of the 7* — > V impact factor and the BFKL Green's function in the MS scheme. 
In addition, a procedure to get results independent from the energy and renormalization scales 
has been investigated within NLA approximation. A shortcoming of those approaches is that 
they consider only the linear regime of the QCD dynamics and nonlinear effects associated 
to the saturation physics are disregarded. However, double meson production in two photon 
interactions at high energies offers an ideal opportunity for studying the transition between 
the linear and saturation regimes since virtualities of both photons in the initial state can 
vary as well as the vector mesons in the final state. In the interaction of two highly virtual 
photons and/or double heavy vector meson production we expect the dominance of hard physics 
(linear regime). On the opposite case, characterized by double light vector meson production 
on real photons scattering, the soft physics is expected to be dominant. Consequently, for an 
intermediate scenario we may expect that the main contribution comes from semi-hard physics, 
determined by saturation effects. 

In this paper we derive the main formulas to describe the double meson production in the 
dipole picture and analyze the double meson production in two photon interactions. We consider 
three cases of physical and phenomenological interest: (a) the interaction of real photons, and 
the interaction of virtual photons with (b) equal and (c) different virtualities. In all cases 
we calculate the forward differential cross section for the pp, pJ/^ and J/tyJ/ty production. 
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Moreover, we present a comparison between the linear and nonlinear predictions and estimate 
the contribution for distinct photon polarizations. 



2 Basic Formulas 



2.1 Double meson production in the dipole picture 

Let us introduce the main formulas concerning the vector meson production in the color dipole 
picture. First, we consider the scattering process 77 — > V\ V 2 , where Vi stands for both light 
and heavy mesons. At high energies, the scattering process can be seen as a succession on 
time of three factorizable subprocesses: i) the photon fluctuates in quark- ant iquark pairs (the 
dipoles), ii) these color dipoles interact and, iii) the pairs convert into the vector mesons final 
states. Using as kinematic variables the 7*7* c.m.s. energy squared s = W 2 = (p + q) 2 , where 
p and q are the photon momenta, the photon virtualities squared are given by Q\ = —q 2 and 
Q\ = —p 2 . The X12 variable is defined by 

T Qi + Ql + M 2 Vi +M 2 2 

W 2 + Q 2 + Q 2 ■ W 

The corresponding imaginary part of the amplitude at zero momentum transfer reads as 

h,h n ' n 

x J dz 2 d 2 r 2 *l fi (z 2 , r 2 , Ql) ^%{z 2 , r 2 ) a dd {x 12 , 77, r 2 ) , (2) 



where \1/ 7 and ty Vi are the light-cone wavefunctions of the photon and vector meson, respectively. 
The quark and antiquark helicities are labeled by h, h, n and n and reference to the meson 
and photon helicities are implicitly understood. The variable 77 defines the relative transverse 
separation of the pair (dipole) and z 1 (1 — z±) is the longitudinal momentum fractions of the 
quark (antiquark). Similar definitions are valid for r 2 and z 2 . The basic blocks are the photon 
wavefunction, \1/ 7 , the meson wavefunction, L , and the dipole-dipole cross section, add- 
In the dipole formalism, the light-cone wavefunctions ^ h ^(z, r) in the mixed representation 
(r, z) are obtained through two dimensional Fourier transform of the momentum space light- 
cone wavefunctions ^/^^(z, k) (see more details, e.g. in Refs. [21, 22, 23]). The normalized 
light-cone wavefunctions for longitudinally (L) and transversely (T) polarized photons are given 
by: 

= ifv^e/Ml-^^, (3) 



K { r ±] ^ r ) = Ml^ee f ie^(z5 h± ^-(l-z)5 h ^ h ±)d r + m f 5 h± - h± , (4) 



2vr 



Kn(er) 



2tt 
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where e 2 = z(l — z)Q 2 



m 



f 



The quark mass rrif plays a role of a regulator when the 



photoproduction regime is reached. Namely, it prevents non-zero argument for the modified 
Bessel functions K 0tl (er) towards Q 2 — > 0. The electric charge of the quark of flavor / is given 
by e e/. 

For vector mesons, the light-cone wavefunctions are not known in a systematic way and 
should be modeled. The simplest approach assumes a same vector current as in the photon case, 
but introducing an additional vertex factor. Moreover, in general the same functional form is 
chosen for the scalar part of the meson light-cone wavefunction. Here, we follows the analytically 
simple DGKP approach [24]. In this particular approach, one assumes that the dependencies 
on r and z of the wavefunction are factorised, with a Gaussian dependence on r. Its main 
shortcoming is that it breaks the rotational invariance between transverse and longitudinally 
polarized vector mesons [25]. However, as it describes reasonably the HERA data for vector 
meson production, as pointed out in Refs. [23, 26], we will use it in our phenomenological 
analysis. The DGKP longitudinal and transverse meson light-cone wavefunctions are given by 
[24], 



*(! - z ) K-h cfntr- fd z ) ex P 



— uj 2 l r 2 



(5) 



loo^ re 



±i6 



m v 



[z&h±,h* - (! - z ) S h T ,h±\ + 



m f 



m v 



h±,h± 



X /oxr- friz) exp 



—ijj\r 2 



(6) 



where e/ is the effective charge arising from the sum over quark flavors in the meson of mass 
my. The following values if = l/y/2 and 2/3 stand for the p and J/ty mesons, respectively. 
The coupling of the meson to electromagnetic current is labeled by f v = 3 m v T e + e - /Ana 2 m 
(see Table 1). The function fr,L(z) is given by the Bauer-Stech-Wirbel model [27]: 



fr,L{z) = N t ,l \jz{l - z) exp 



-m 2 v (z - l/2f 
2uj< l 



(7) 



The meson wavefunctions are constrained by the normalization condition, which contains the 
hypothesis that the meson is composed only of quark-antiquark pairs, and by the electronic 
decay width T v ^ e + e -. Both conditions are respectively given by [28, 22], 



(8) 



h,h 



E 

Kh 



j^y 2 J{T ^ 7) [z(l ~ Z )Q 2 + k 2 + m)\ *Xfi(k, z)*l*- h (k, z) = efymy (e; ■ e v ) . (9) 



The constraints above, when used on the DGKP wavefunction, imply the following relations 
[23], 
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(10) 
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V(m v ) 


ev 


fv 


UJt 


Mr 




M L 


MeV 




[GeV] 


[GcV] 




[GeV] 




p(770) 


l/y/2 


0.153 


0.218 


8.682 


0.331 


15.091 


J/V (3097) 


2/3 


0.270 


0.546 


7.665 


0.680 


19.350 



Table 1: Parameters and normalization of the DGKP vector meson light-cone wavef unctions. 
Results obtained using quark mass values from the saturation model (see text). 

f 1 , ^ , s f 1 2 \z 2 + (1 - z) 2 ] ujl + ml , N 



where 



/ L = ^ dzzH\-zffl(z), (12) 

° lfc [^ + (l- 2 )>- +m?/IM (13) 

my 



The relations in Eq. (10) come from the normalization condition, whereas the relations in Eq. 
(11) are a consequence of the leptonic decay width constraints. The parameters u>t,l and Mt,l 
are determined by solving (10) and (11) simultaneously. In Table 1 we quote the results which 
will be used in our further analysis. To be consistent with the saturation models, which we 
will discuss further, we have used the quark masses m U:diS = 0.14 GeV and m c = 1.5 GeV. We 
quote Refs. [23, 26] for more details in the present approach and its comparison with data for 
both photo and electroproduction of vector mesons. 

Finally, the imaginary part of the forward amplitude can be obtained by putting the expres- 
sions for photon and vector meson (DGKP) wavef unctions, Eqs. (3-4) and (5-6), into Eq. (2). 
Moreover, summation over the quark/antiquark helicities and an average over the transverse 
polarization states of the photon should be taken into account. In order to obtain the total 
cross section, we assume an exponential parameterization for the small \t\ behavior of the am- 
plitude. After integration over \t\, the total cross section for double vector meson production 
by real/virtual photons reads as, 

a (77 -> Vi V 2 ) ~ 



By 1 y 2 dt 



[lmA(s, t = 0)] 

16"7T By l V 2 



(14) 



where (5 is the ratio of real to imaginary part of the amplitude and By l y 2 is the slope parameter. 



2.2 Dipole-dipole cross section in the saturation model 

The dipole formulation has been extensively used in the description of inclusive and diffractive 
processes at HERA in an unified way. The basic quantity is the dipole-proton cross section a^p, 
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which contains all information about the target and the strong interaction physics. In general, 
the saturation models [5, 6, 7, 8] interpolate between the small and large dipole configurations, 
providing color transparency behavior, a div ~ r 2 , at r l/Qsat an d constant behavior at 
large dipole separations r > 1/Q sa t- The physical scale which characterizes the transition 
between the dilute and saturated system is denoted saturation scale, Q 2 at oc x~ x , which is 
energy dependent. Along these lines, the phenomeno logical saturation model proposed by 
Golec-Biernat and Wusthoff (GBW) [5] resembles the main features of the Glauber-Mueller 
resummation. Namely, the dipole cross section in the GBW model takes the eikonal-like form, 



Its phenomenological application has been successful in a wide class of processes with a photon 
probe. Although the GBW model describes reasonably well the HERA data, its functional form 
is only an approximation of the theoretical nonlinear QCD approaches [3, 4]. The parameters 
of the model are <7 = 29.12 mb, light quark masses (u,d,s) rrif = 0.14 GeV, charm mass 
m c = 1.5GeV. Moreover, the saturation scale is given by Qsat = (xo/x) x ^ 2 , with the parameters 
x — 4.1 • 10" 5 and A = 0.277. In Ref. [8] a parameterization for the dipole cross section was 
constructed to smoothly interpolate between the limiting behaviors analytically under control: 
the solution of the BFKL equation for small dipole sizes, r <C 1/Q s &t{ x ), and the Levin- Tuchin 
law [29] for larger ones, r ^> 1/Q sa ,t(x). A fit to the structure function F 2 (x, Q 2 ) was performed 
in the kinematical range of interest, showing that it is not very sensitive to the details of the 
interpolation. The dipole cross section was parameterized as follows, 



where the expression for rQ sat (x) > 2 (saturation region) has the correct functional form, as 
obtained either by solving the Balitsky-Kovchegov (BK) equation [3], or from the theory of 
the Color Glass Condensate (CGC) [4]. Hereafter, we label the model above by IIM. The 
coefficients a = 1.8 and b = 1.13 are determined from the continuity conditions of the dipole 
cross section at rQ SSLt (x) = 2. The coefficients 7 sat = 0.63 and k = 9.9 are fixed from their 
LO BFKL values. In our further calculations it will be used the parameters R p = 0.641 fm, 
A = 0.253, x = 0.267 x 10~ 4 and Af = 0.7, which give the best fit result. It is important to 
emphasize that the GBW and IIM saturation models are suitable in the region below x = 0.01 
and the large x limit needs still a consistent treatment. At ep collisions the dipole-proton cross 
sections should be supplemented by a threshold factor (1 — x) TlthICS , with n t hrcs = 5, which is 
directly associated with the number of spectators at x ~ 1 (n t hrcs = 2n spcct — 1). 

Following Ref. [9] we can extend the saturation model, originally proposed to describe 
ep collisions, to two-photon interactions at high energies. The basic idea is that the dipole- 
dipole cross section <Jdd( x i2, r i, r 2) has the same functional form as the dipole-proton one and 
is expressed in terms of an effective radius r cS , which depends on ri and/or r 2 . Consequently, 
we have that [9], 




(15) 




A/"o (TC^p— — ; , forrQ sat (:r)<2, 
1 - exp -a In 2 (brQ sat ) , for rQ sat (x) > 2 , 



(16) 




(17) 
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Figure 1: Dipole-dipole cross sections: (a) Comparison between the predictions for GBW and 
IIM models at different values of x\2- (b) Comparison between their resummed predictions and 
the corresponding linear limits. 



and 

ln(2/r cff Q sat ) 



off 1 («»,r„) = »A M " , ' ■ fOT »*««<*») S 2 . (18) 

[ 1 - exp -a In 2 (6r cff Q sat ) , for r eS Q sat (x 12 ) > 2 , 

where the x\2 variable is given by the Eq. 1 and &q = |<7o, with <7o the same as in Refs. [5] and 
[8] and referred above. The last relation can be justified in terms of the quark counting rule. 
In the two-photon case, the resulting quark masses are slightly different as it was found in Ref. 
[9]: m/ = 0.21 GeV for light quarks and m c = 1.3 GeV for charm. In Ref. [9] three different 
scenarios for r e g has been considered, with the dipole-dipole cross section presenting in all cases 
the color transparency property (cr dd — > for r\ — > or r 2 — > 0) and saturation (cr^ — > <r ) for 
large size dipoles. We quote also Ref. [30]) for interesting discussions on the effective radius 
and its consequences in hadron-hadron interactions. In what follows, we use the model I from 
[9], where r 2 ff = r\r\j{r\ + r 2 ), which is favored by the 7*7* and data. We have tested 
the sensitivity of the result to a different prescription, r 2 ff = min (r 2 , r 2 ,) (named model II in 
Ref. [9]). Its deviation from model I is quite large for p production and almost insensitive for 
the mixed pJ/^ production. For double J/ty production the deviation is considerably larger 
than the mixed one. However, the difference concerns only to the overall normalization and no 
change is seen in the energy behavior. Moreover, in order to extend the dipole model to large 
X12 it is necessary to take into account threshold correction factors which constrain that the 
cross section vanish when x 12 — > 1 as a power of 1 — X\i. As in Ref. [9], we multiply the dipole- 
dipole cross section by the factor (1 — X12) 5 . A comment is in order here. One shortcoming 
of the GBW model is that it does not contain the correct DGLAP limit at large virtualities. 
Consequently, we may expect that its predictions are only valid at small values of the photon 
virtualities. Therefore, in what follows we only consider photon virtualities up to 10 GeV 2 . 
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(a) (b) 

Figure 2: A comparison between the predictions for (a) total cross section and (b) charm 
production in two real photon collisions considering the GBW and IIM dipole-dipole cross 
sections. The QPM contribution has been added. Data from OPAL and L3 Collaborations. 



In Fig. 1 (a) we present the dependence of the two dipole-dipole cross sections, Eqs. (17) and 
(18), as a function the effective radius r eS at different values of a; 12 (a; 12 = lCT n , n— 1,2, 3, 4). 
We have that at small values of r e g their predictions are similar, while they differ approximately 
15 % at large r e fj and small values of x\2- In order to emphasize the importance of the saturation 
effects, in Fig. 1-b we present a comparison between the full predictions of the GBW and IIM 
dipole-dipole cross sections and their linear limits. We have denoted by resummed the curves 
with the complete expressions in Eqs. (17-18) and by linear their approximations in the limit 
of small dipoles. Namely, for the linear case one has cx2f w oc <3"o(<3sat r eff/4) and for the IIM 
model we just take the extrapolation of Eq. (18) for r c g < 2/Q SSut . We have that at r e g rs 0.2 
fm the linear and resummed predictions from the GBW model start to be different. On the 
other hand, in the IIM case, this difference starts at r g rs 0.5 fm. Consequently, the transition 
between the linear and saturation regimes is distinct in the GBW and IIM models. 

In next section we will compare the predictions for the double meson production coming 
from different models for the dipole-dipole cross section. However, the extension of the IIM 
model for photon-photon interactions has not been considered before. For sake of completeness, 
we compare the predictions of the GBW and IIM models for the specific cases of the total 77 
cross section and heavy quark production. The analysis of the heavy quark production is 
motivated by its strict relation with the double J/ty production. The GBW model has already 
been considered in Ref. [9], while the IIM analysis is the first one in literature. In Fig. 2 (a) we 
present a comparison between the predictions of the GBW and IIM models for the total cross 
section and the OPAL and L3 experimental data [31, 32]. Following [9] we also include the 
QPM and reggeon contribution and assume the model I for the effective radius. We have that 
the GBW and IIM predictions are similar, describing the current experimental data quite well. 
It should be noticed that parameters for IIM have not been adjusted in order to fit two-photon 
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data as done for GBW. Furthermore, we compute the cross section for charm production in 
the reaction 77 — > cc, considering real photons. The results are presented in Fig. 2 (b) for two 
prescriptions of the effective radius (model I and II referred before) and compared with L3 data. 
The low energy quark box contribution (QPM) has been added. An additional contribution, 
which we do not include, is the resolved (single and double) piece to the charm cross section, 
which reaches 30 % of the main contribution at high energies. As already verified in Ref. [9], 
both prescriptions for the effective radius provide reasonable description of the data when GBW 
model is considered. On the other hand, in the IIM model, the prescription I for the effective 
radius gives a better description, with the model II overestimating the L3 data at high energies. 
Moreover, the IIM model implies a stronger energy dependence of the heavy quark production 
cross section than the GBW prediction. This behavior should also be present in other processes 
characterized by a hard scale as, for instance, the interaction of two highly virtual photons or 
double heavy vector meson production. In what follows we only will consider the model I for 
the effective radius. 



3 Results 

In order to calculate the total cross section for double vector meson production given in Eq. 
(14) it is necessary to specify the value of the slope parameter By x y 2 . As this quantity is 
not well constrained, in what follows we only will present our predictions for the energy and 
virtuality dependence of the forward differential cross section dcr( - 77 7^ yi v ^ _ ■ This should 

be enough for the present level of accuracy. We start by the study of the scattering of two real 
photons, investigating its dependence on energy and on the mesons mass. After, we consider 
the scattering of virtual photons and investigate the symmetric (Q\ = Q%) and asymmetric 
{Q\ oc aQ\ with a >> 1) cases. In addition, we estimate the magnitude of the contribution 
of the distinct polarizations for the total cross sections. Finally, we discuss the size of parton 
saturation effects in the production of different mesons. 



3.1 Double meson production on real photons interactions 

Let us start our analyzes considering the double meson production in two real photon scattering. 
In Fig. 3 we present the forward differential cross sections for the representative cases of double 
p, pJ/^ and double J/\I/ production in the energy range 50 GeV < W 77 < 10 3 GeV. The curves 
are presented for the two models of the dipole-dipole cross section given in Eqs. (17) and (18). 
Bold curves stand for GBW model and thin curves for IIM model. The forward differential 
cross section is sizeable in the double p case, being of order 20-40 nb/GeV 2 in the energy range 
considered. The mixed pJ/^ production is the second higher rate reaching 4 — 40 pb/GeV 2 , 
whereas double J/^f production is quite low. The deviations between the GBW and IIM 
models are large for the double p production, with IIM results being a factor 10 below GBW 
one at ~ 1 TeV. The origin of this discrepancy is not clear, since there is no evidence 
for strong deviations in the 77 case shown in Fig. 2-a. Probably, deviations could come from 
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Figure 3: Energy dependence of the forward differential cross section for double vector meson 
production considering real photons interactions {Q\ = Q\ = GeV 2 ). Bold (thin) curves are 
the results for GBW (IIM) parameterization for dipole-dipole cross section. 



the different weights given by the wavefunctions in each case. This subject requires further 
investigations. For the double J/\l/ production, the IIM prediction overestimates the GBW 
by a factor 4, which agrees with the expectation which comes from our previous results for 
heavy quark production (See Fig. 2-b). On the other hand, in the pJ/^ case, the results 
are equivalent at low energies but differ by a factor two at 1 TeV, with the GBW prediction 
being greater than the IIM. These features can be qualitatively understood in terms of the 
scales involved in the process. As we discussed before, the IIM dipole-dipole cross section has a 
relatively faster transition to saturation in comparison with GBW and underestimates it by a 
factor of 20-30% at small-Xi2- In double p and mixed vector meson production the typical scale 
is given by the light meson mass p 2 = 2M 2 or the sum of light-heavy meson p 2 = (M 2 + M 2 ^). 
Therefore, the double p process is dominated by a relatively soft scale and saturation effects 
should be important, whereas the mixed production is characterized by a semihard scale which 
is still sensitive to saturation effects. On the other hand, in the double J/ty production the 
typical scale is sufficiently hard, p 2 = 2M 2 ^. Therefore, we expect a larger contribution of 
small dipoles leading to a cross section with higher magnitude. 

In order to analyze the energy dependence of the forward differential cross section we have 
performed a simple power-like fit in the energy interval 50 < W 77 < 10 3 GeV in the form 
_ oc W« r For the double p production one obtains a = 0.4 (0.08) for GBW (IIM) 
model. In Regge phenomenology, this corresponds to an effective Pomeron intercept of order 
a c p ~ a/4 = 0.1 (0.02) for GBW (IIM) parameterizations, which is clearly a soft behavior. 
This fact shows that the IIM model contains stronger saturation effects in contrast to GBW 
one in the case of dominantly soft scales. In the mixed production, the effective power increases 
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Figure 4: Energy dependence for double vector meson production considering virtual photons 
considering: (a) equal virtualities {Q\ = Q\ = 10 GeV 2 ) and (b) distinct virtualities (Q\ = 
and Ql = 10 GeV 2 ). 



to a = 0.96 (0.65) for GBW (IIM) model and the difference is not too sizeable as in the p case. 
For the double J/^/ production, a = 1.06 (0.9), which implies a e p ~ 0.27. Therefore, one has 
a hard Pomeron behavior in the case where the heavy meson mass (p 2 = 2Mj,^,) is present in 
the problem. Thus, as expected from the phenomenology of ep collisions, the saturation model 
for double vector meson production in 77 interactions is able to consistently connect the soft 
behavior when a non-perturbative scale is involved and the hard Pomeron expectations when 
a perturbative scale is present. 

Let us now compare our results with those obtained in other approaches [12, 13, 14, 15, 
16, 17]. Initially, let us consider the previous calculations within the color dipole picture [12]. 
In Ref. [12] there are estimations of the total cross section for double meson production. 
Our predictions underestimate those results by a factor ten for double p and a factor one 
hundred for the other mesons. In this comparison we have used the following values for the 
slope parameters: B pp = 10 GeV~ 2 , B pJ /q, = 5 GeV -2 and Bj/^j/^ = 0.44 GeV -2 , which 
are taken from our recent investigations on double meson production in Refs. [15, 16]. These 
deviations are probably due to the different dipole-dipole cross section, distinct choices for the 
quark masses and uncertainties in the determination of the slope parameter. For instance, the 
dipole-dipole cross section in Ref. [12] behaves as add oc r \ r \ for small dipoles and add oc r\r\ 
for large dipoles, which overestimate the integration on dipole sizes in comparison with the 
dipole-dipole cross sections presented here. Namely, one has add ^ r esQsat f° r dipoles having 
transverse size r eff < 1/Q sat and add <5"o for dipoles of size r eS > l/Qsat- Furthermore, the 
pJ/'ty production in 77 processes has also been estimated in Refs. [14, 16]. There, the differential 
cross section was estimated in a similar way as the elastic J/\I/ photoproduction off the proton 
[33]. Our results agree with these predictions, with a behavior similar to those obtained using 
the GRS(LO) parameterization for the gluon distribution on the p meson. This process was 
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Figure 5: Energy dependence for double p production at equal photon virtualities, Q 2 = 1 and 
Q 2 = 10 GeV 2 . The contributions of the different polarizations (TT, LL, TL/LT) are explicitly 
presented. See discussion in text. 



also estimated in Ref. [17] using the non-forward solution of the BFKL equation. Our results 
are smaller than estimations obtained in Ref. [17]. This is expected since the saturation effects 
modifies significantly the cross section of this semi-hard process. Moreover, our results for the 
double p production agree with those obtained in Ref. [16] assuming the pomeron-exchange 
factorization. However, as its predictions are strongly dependent on the assumptions present 
in the calculations of the double J/^f and pJ/^ production (See Table 1 in Ref. [16]), a 
direct comparison is not very illuminating. This process also was analyzed in Ref. [17], but a 
direct comparison is not possible because only the hard contribution (|£| > 1 GeV 2 ) has been 
estimated. 



3.2 Double meson production on virtual photons interactions 

Let us now consider the double vector meson production when we have the interaction of 
virtual photons. In Fig. 4 (a) we present the predictions of the GBW model for the energy 
dependence considering that the incident photons have equal virtualities (Q 2 = 10 GeV 2 ). We 
have that the forward differential cross section decreases when the virtuality and/or the total 
mass of the final state is increased. The differential cross sections present a behavior similar 
on energy, independently of the meson mass. This is due the sufficiently hard scale for these 
processes given by p 2 = 2Q 2 + M v + M V2 , which is basically determined by the high photon 
virtuality, since 2Q 2 > My. This also explains the proximity between the pJ/^ and double 
J/^> predictions, in contrast with obtained for the real photon interactions. A power-like fit 



to the differential cross section in the form 



d<? y 1 v 2 
dt 



oc gives a = 1.01, 1.08, 1.1 for 
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double p, pJ/^ and double J/^f, respectively. Our result for double p is consistent with the 
NLA BFKL calculation using BLM scale fixing presented in Ref. [19]. 

In Fig. 4 (b) we present our predictions for double meson production considering unequal 
photon virtualities. We consider the limit case of real photon scattering on a deeply virtual 
partner, namely Q\ = and Q\ = 10 GeV 2 . Now, the typical scale is given by p 2 = Q 2 + 
M Vi + My 2 . In our computation of the mixed production we take the following statement for 
the photon virtualities: Q\ corresponds to the photon transforming into p and Q\ corresponds 
to the photon transforming into J . Notice that the final cross section should be given by 
f 7[p(g 1 )J/*(g 2 )orp(Q 2 )J/*(Q 1 )] = <7b(Qi)J/tt(Q 2 )] + ( j[p(Q 2 )J/*(Q 1 )]. We have that the 
behavior of the different predictions are similar those obtained in Fig. 4 (a), with the energy 
dependence for pJ/^ and double J/\I/ production being almost identical those obtained in the 
symmetric case. The main difference occurs for double p production, which has its energy 
dependence strongly modified by saturation effects due to the small value of p 2 present in the 
problem. A power-like fit to the differential cross section in the form duy ^ 2 t n oc W" 7 gives 
a = 0.65, 1.03, 1.08 for double p, pJ/^ and double J/^f, respectively. 



=o 



Using the dipole approach the contribution of the different polarizations for the forward 
differential cross section can be directly estimated. Let us start considering the double p 
production at equal virtualities of photons. We consider the following illustrative cases, Q\ = 
Q\ — 1 GeV 2 and Q\ = Q\ = 10 GeV 2 . These choices allow to observe the dependence of 
each contribution on virtuality. The results are shown in Fig. 5. The transverse piece (TT) 
is labeled by dotted curves, longitudinal piece by dashed curves, mixed transverse-longitudinal 
(TL or LT) by long dashed curves and the total cross section (summation over polarizations) 
by the solid curves. In case of production of a same vector meson, the TL and LT pieces 
contribute equally, TL = LT. For virtualities Q 2 = 1 GeV 2 , the transverse content dominates, 
followed by the LT/LT and LL pieces. Longitudinal content is a quite small contribution, which 
is consistent with the longitudinal wavefunction to be proportional to photon virtuality, which 
vanishes when Q 2 — > 0. A completely different situation occurs when the virtualities increase 
to Q 2 = 10 GeV 2 . In this case the longitudinal piece is dominant, followed by the LT/LT and 
transverse parts. This is consistent with the ratio o"l/o" T > being Q 2 -dependent in the light 
meson photoproduction (See e.g. Ref. [25]). 

A similar analysis can be made for the double J/^f production (see Fig. 6). We take 
the same virtualities for the virtual photons and equal notation as before. For virtualities 
Qi = Q\ = 1 GeV 2 , the transverse content dominates, followed by the LT/LT and LL pieces. 
As in the double p case, the longitudinal contribution is quite small. The total contribution 
is determined completely by the transverse contribution, with other pieces being negligible. A 
completely different situation occurs when virtualities increase to 10 GeV 2 in contrast with the 
p case. The pattern remains the same as for Q 2 = 1 GeV 2 , with transverse piece still dominant, 
followed by LT/TL and LL pieces. The total contribution is slightly larger than the transverse 
one. 

Finally, let us investigate the dependence on virtuality at fixed energy of the forward dif- 
ferential cross section. We take the representative energy of = 500 GeV. In Fig. 7 we 
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Figure 6: Energy dependence for double J/ty production at equal photon virtualities, Q 2 = 1 
and Q 2 = 10 GeV 2 . The contributions of the different polarizations (TT, LL, TL/LT) are 
explicitly presented. See discussion in text. 



present the dependence of the forward differential cross section on the ratio R = Q\jQ\ at 
fixed Q\. We consider the typical values Q\ = 1 and Q\ = 10 GeV 2 . In all cases, the cross 
section decreases as Q\ increases, presenting finite values towards R = 0. It should mentioned 
that this interpolation can not be obtained in the BFKL approach in view of the lack of a 
scale in the process. In our case, the saturation scale provides the semihard scale. In order 
to investigate the quantitative behavior on R at intermediate virtualities, we adjust the curves 
with the simple exponential parameterization for R > 1 in the form day ^ 2 oc e~^ R . This 

procedure gives (3 = 4.4, 2.6, 2.9 for double p, double J/^f and pJ/^, respectively. The results 
follow the typical saddle-point BFKL solution in the region of R > 5 GeV 2 , namely the cross 



noticed that our definition for R is slightly different from that reference 



sections behave as — 4- 1 



oc r&rfi exp(— (3 In 2 R) as computed in Ref. [19]. It should be 



3.3 Investigating saturation effects 

Let us now investigate the magnitude of saturation effects in the differential cross section 
comparing the results using the small r e fj approximation for the dipole-dipole cross section 
with the complete expression including the transition for the saturation regime (See discussion 
in Sect. 2.2). We restrict our analysis to a comparison between the linear and saturation model 
predictions for double p and J/^f production. The results are shown in Fig. 8, where solid lines 
stand for the resummed calculations and dot-dashed one for the linear approximation. It should 
be noticed that the saturation scale is different for each meson because Q 2 at ~ (20/212) ' 3 anc ^ x vi 
depends on the meson mass as defined in Eq. (1). For the most striking case, in the production 
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Figure 7: Dependence of the forward differential cross section on the ratio R = Q\jQ\ of photon 
virtualities at Q\ fixed. 

by two real photons, the saturation scale for p reaches Q 2 at ~ 2.5 GeV 2 whereas stays as Q 2 at ~ 1 
GeV 2 for J/ty. Therefore, the saturation scale is higher for p than for J/^> up to intermediate 
virtualities. Let us start discussing double J/^f production [see Fig 8 (b)], which is one typical 
hard process characterized by a hard scale given by p 2 {Q 1 M v ) = 2 (Q 2 + M 2 j^). Consequently, 
we may expect that a perturbative description to be valid even in the real photon limit, Q 2 —>■ 0, 
and that the contribution from the saturation effects to be small while p 2 3> Q 2 at . However, 
as the saturation scale grows with the energy, the saturation effects become important at large 
energies. This is the reason we observe a difference between linear and resummed predictions at 
W ~ 1 TeV in the real photon case. These results indicate that double J/\P production is not 
strongly modified by saturation corrections for energies smaller than 1 TeV. On the other hand, 
in the double p production the situation changes drastically. For real photons, it is a typical 
soft process and, therefore, in this case the linear and saturation predictions are very distinct. 
Now, the scale is given by p 2 (Q, My) = 2 (Q 2 + M 2 ) , which should be treated carefully due to 
the small meson mass. The results are presented in Fig. 8 (a). In the real photon scattering, 
one has p 2 = 2M 2 m 1 GeV 2 and therefore saturation effects are increasingly important as 
p 2 < Qsat(^ry)> which explain the reason for the cross section to be reduced by two orders of 
magnitude at W 77 ~ 1 TeV. At Q 2 = 10 GeV 2 , the scenario is different since p 2 = 2Q 2 > Q 2 at . 
Therefore, the resummed prediction is similar to the linear but corrections for large energies 
are still important. In view of discussions above, double p production becomes an ideal place 
to probe the saturation physics. 
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Figure 8: Comparison between the linear and saturation predictions for the energy dependence 
of the forward differential cross section: (left panel) double p production and (right panel) 
double J/ty production. 

4 Summary 

In this paper we have extended the dipole picture for the double vector meson production, 
l*(Qi)l*(Q2) —> V1V2, and calculated the forward differential cross section assuming that the 
dipole-dipole cross section can be modeled by a saturation model. We have analyzed the 
energy and virtuality dependence and investigated the magnitude of saturation effects. It is 
found that the effective power on energy is directly dependent on the typical momentum scale 
for the process, p? = Q\ + Q\ + M v + My 2 , which is different for distinct meson pair and 
photon virtualities. Saturation effects are important for double p production on real photons, 
whereas is small for processes containing J/ty and/or large photon virtualities. It is shown 
the contribution of the distinct polarizations and their regions of dominance for each meson 
pair. The results are consistent with expectations from electroproduction of vector mesons. 
The dependence on virtuality has been investigated using the analysis on the ratio R = Q\jQ\- 
The results are qualitatively in agreement with previous predictions obtained using dipole or 
NLO BFKL approaches. Our results demonstrate that double meson production in two photon 
interactions at high energies offer an ideal opportunity for the study of the transition between 
the linear and saturation regimes. 



16 



Acknowledgments 



VPG would like to thanks W. K. Sauter for informative and helpful discussions. One of us 
(M. Machado) thanks the support of the High Energy Physics Phenomenology Group, GFPAE 
IF-UFRGS, Brazil. The authors are grateful to C. Marquet for his valuable comments and 
suggestions. This work was partially financed by the Brazilian funding agencies CNPq and 
FAPERGS. 



References 

[1] L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); E. A. Kuraev, L. N. Lipatov, V. S. 
Fadin, JETP 45, 1999 (1977); I. I. Balitskii, L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 
(1978). 

[2] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15, 438 (1972); G. Altarelli and G. 
Parisi, Nucl. Phys. B126, 298 (1977); Yu.L. Dokshitzer, Sov. Phys. JETP 46, 641 (1977). 

[3] I. Balitsky, Nucl. Phys. B 463, 99 (1996); Y. V. Kovchegov, Phys. Rev. D 60, 034008 
(1999); Phys. Rev. D 61, 074018 (2000). 

[4] L. D. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 (1994); E. Iancu, A. Leonidov, 
L. McLerran, Nucl. Phys. A 692, 583 (2001); E. Ferreiro, E. Iancu, A. Leonidov, L. 
McLerran, Nucl. Phys. A 703, 489 (2002); J. Jalilian-Marian, A. Kovner, L. McLerran 
and H. Weigert, Phys. Rev. D 55, 5414 (1997); J. Jalilian-Marian, A. Kovner and H. 
Weigert, Phys. Rev. D 59, 014014 (1999), ibid. 59, 014015 (1999), ibid. 59 034007 (1999); 
A. Kovner, J. Guilherme Milhano and H. Weigert, Phys. Rev. D 62, 114005 (2000); H. 
Weigert, Nucl. Phys. A703, 823 (2002). 

[5] K. Golec-Biernat and M. Wusthoff, Phys. Rev. D 60, 114023 (1999); Phys. Rev. D 59, 
014017 (1998). 

[6] J. Bartels, K. Golec-Biernat and H. Kowalski, Phys. Rev D66 (2002) 014001. 

[7] H. Kowalski and D. Teaney, Phys. Rev. D 68, 114005 (2003). 

[8] E. Iancu, K. Itakura and S. Munier, Phys. Lett. B 590, 199 (2004). 

[9] N. Timneanu, J. Kwiecinski, L. Motyka, Eur. Phys. J. C 23, 513 (2002). 

[10] N. N. Nikolaev and B. G. Zakharov, Z. Phys. C49, 607 (1991); Z. Phys. C53, 331 (1992); 
A. H. Mueller, Nucl. Phys. B415, 373 (1994); A. H. Mueller and B. Patel, Nucl. Phys. 
B425, 471 (1994). 

[11] N. N. Nikolaev, J. Speth and V. R. Zoller, Eur. Phys. J. C 22, 637 (2002); J. Exp. Theor. 
Phys. 93, 957 (2001) [Zh. Eksp. Teor. Fiz. 93, 1104 (2001)]. 



17 



[12] A. Donnachie, H. G. Dosch and M. Rueter, Phys. Rev. D 59, 074011 (1999) 

[13] J. Kwiecinski, L. Motyka, Phys. Lett. B 438, 203 (1998); Acta Phys. Pol. B 30, 1817 
(1999). 

[14] L. Motyka, B. Ziaja, Eur. Phys. J. C 19, 709 (2001). 

[15] V. P. Goncalves and M. V. T. Machado, Eur. Phys. J. C 28, 71 (2003) 

[16] V. P. Goncalves and M. V. T. Machado, Eur. Phys. J. C 29, 271 (2003). 

[17] V. P. Goncalves and W. K. Sauter, Eur. Phys. J. C 44, 515 (2005). 

[18] B. Pire, L. Szymanowski and S. Wallon, Eur. Phys. J. C 44, 545 (2005) 

[19] R. Enberg, B. Pire, L. Szymanowski and S. Wallon, Eur. Phys. J. C 45, 759 (2006) 

[20] D. Y. Ivanov and A. Papa, Nucl. Phys. B 732, 183 (2006) 

[21] V. Barone and E. Predazzi, High-Energy Particle Diffraction, Springer- Verlag, Berlin Hei- 
delberg, (2002). 

[22] S. Munier, A. M. Stasto and A. H. Mueller, Nucl. Phys. B 603, 427 (2001). 

[23] J. R. Forshaw, R. Sandapen and G. Shaw, Phys. Rev. D 69, 094013 (2004). 

[24] H. G. Dosch, T. Gousset, G. Kulzinger and H. J. Pirner, Phys. Rev. D55, 2602 (1997). 

[25] N. N. Nikolaev, Comments Nucl. Part. Phys. 21, 41 (1992); I. P. Ivanov, N. N. Nikolaev 
and A. A. Savin, arXiv:hep-ph/0501034. 

[26] V. P. Goncalves and M. V. T. Machado, Eur. Phys. J. C 38, 319 (2004). 

[27] M. Wirbel, B. Stech and M. Bauer, Z. Phys. C29, 637 (1985). 

[28] S. J. Brodsky and G. P. Lepage, Phys. Rev. D22, 2157 (1980). 

[29] E. Levin and K. Tuchin, Nucl. Phys. B 573, 833 (2000). 

[30] C. Marquet and R. Peschanski, Phys. Lett. B 587, 201 (2004). 

[31] G. Abbiendi et al. [OPAL Collaboration], Eur. Phys. J. C14, 199 (2000). 

[32] M. Acciarri et al. [L3 Collaboration], Phys. Lett. B519, 33 (2001). 

[33] M. G. Ryskin, Z. Phys. C 57, 89 (1993). 



18 



